Pulsars With Jets May Harbor Dynamically Important Accretion Disks 
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(N ■ ABSTRACT 
O 

Q ' For many astrophysical sources with jets, there is evidence for the contemporaneous 

^j- ■ presence of disks. In contrast, pulsars such as the Crab and Vela show jets but have 

not yet revealed direct evidence for accretion disks. Here we show that for such pulsars, 
an accretion disk radiating below detectable thresholds may simultaneously account 
for (1) observed deviations in the braking indices from that of the simple dipole, (2) 
observed pulsar timing ages, and (3) possibly even the jet morphology via a disk 
outflow that interacts with the pulsar wind within, collimating and/or redirecting it. 
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52 1 1. Introduction 

C3 ■ 

Recent X-ray observations have confirmed that collimated jets can be found emanating from 
pulsar engines (Weisskopf et al. 2000; Helfand et al. 2001; Gaensler et al. 2002; Slane 2002; 
Pavlov et al. 2003). These jets morphologically resemble those of young stellar objects (Reipurth 
& Bally 2001), planetary nebulae, (Balick &; Frank 2002) active galactic nuclei (Ferrari 1998) and 
microquasars (Mirabel 2001). Jets and disks are likely associated in these systems, even if the 
role of the disk vs. central object is not fully understood with respect to jet production. Here we 
explore the possibility that there may also be disks around jetted pulsars, and that winds from 
these disks might collimate the pulsar wind. The motivation is that relativistic pulsar winds do 
not seem to self-collimate easily (e.g. Lyubarsky &: Eichler 2001; Bogovalov & Khangoulian 2003; 
Komissarov & Lyubarsky 2003): strong outward electric forces compete with the self-collimating 
collimating hoop stresses. Tsinganos & Bogovalov (2002) simulated the generic efficacy of 
relativistic wind collimation by an ambient non-relativistic wind, but did not dynamically relate 
the external wind to the properties of an underlying disk. 

Michel & Dessler (1981) suggested that disks formed from the stellar collapse offer a physical 
similarity between radio and X-ray pulsars. Such disks can also influence pulsar braking indices 
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and timing ages: the former often deviate from those expected if spin-down were due only to dipole 
radiation, and the latter often differ from inferred actual ages. Various intrinsic mechanisms for 
these deviations have been proposed (Macy 1974; Manchester & Taylor 1977; Chanmugam & Sang 
1989; Beskin et al. (1993); Chubarian et al. 2000; Wu et al. 2003), including an external torque 
provided by a fallback accretion disk (e.g. Menou et al. 2001a; Marsden et al. 2001; Alpar et al. 
2001; Qiao et al. 2003; Shi & Xu 2003). This would imply a weaker pulsar magnetic field than 
that inferred from the measured period P and its derivative P. This has been confirmed for one 
neutron star where the magnetic field was measured from cyclotron lines (Bignami et al. 2003). 

Here we calculate the effect of a wind-emitting accretion disk on the pulsar spin evolution 
by coupling the latter to evolution equations for accretion and radius dependent mass loss. We 
constrain the pulsar magnetic field strength and the initial outer accretion rate (or disk mass) for 
which the solutions match the observed periods, timing ages, and braking indices of Crab and 
Vela. The resulting solutions produce sufficient disk-wind momenta to potentially collimate and 
help illuminate the pulsar winds. We describe our wind and accretion parameterization in section 
2. In Sec. 3 we show that the launch mechanism determines the radial dependence of the mass 
outflow rate. In Sec. 4 we calculate the spin evolution of the pulsar and find best fit solutions for 
Crab and Vela. We compare the pulsar and disk-wind momenta in Sec. 5, and conclude in Sec. 6. 



2. The accretion and disk-wind 

Given the accretion rate at the outer disk radius M out , the time-dependent mass loss from a 
disk-wind launched between radii i? ou t and R < R out from the pulsar can be parameterized as 



MdwOR,*) = Mout(t) 



1 



R 



Rout (t) 
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(Values of M ont (t = to) an d the parameter < p < 1 will be discussed below.) The accretion rate at 
the inner edge of the disk where the torque is exerted, is given by Mj n (t) = M out (t) — M^ w (Ri n ,t). 
The disk is allowed to penetrate the magnetosphere up to the magnetospheric radius 1 
i? in ~ 2.55 x 10 8 M^M^B^ 7 cm (£i 2 10 12 G is the pulsar's magnetic field and M NS ,iM its 
mass) where the torque exerted by the magnetic field on the disk is of order the viscous torque. 

To constrain M out (t), we appeal to previous work. Chevalier (1989) studied fallback accretion, 
but without angular momentum. Here we adopt the Menou et al. (2001a; see also Chatterjee et al. 
2000) analogy between fallback accretion and accretion of the torus of gas formed by a disrupted 
star (Cannizzo et al. 1990). Cannizzo et al. (1990) showed that there is a transient spreading 
phase of duration to during which the accretion is nearly constant, M out (t) ~ M ut(^o) but after 



lr The surface where the pulsar field develops a significant toroidal component (facilitating a Poynting flux outflow) 
likely has an hourglass shape, being pinched at the disk midplane. Field lines extending from the pulsar to large 
distances would develop toroidal fields closer to the light cylinder, due to the rapid density drop away from the disk. 
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which the accretion rate declines as a power-law with time, M out (t) = M out (to)[t/to]~ 5 , and the 
initial accretion rate depends linearly on the product of the viscosity parameter a ss (Shakura & 
Sunyaev 1973) and the initial mass of the disk. They also found that 5 = 19/16 provides a good 
fit to the late time evolution and that the outer edge of the disk then evolves as -R ou t oc (t/to) 3 ^ 8 . 
The initial i? ut is determined by the angular momentum of the presupernova star and the amount 
of ejected material. Menou et al. (2001b) estimate -R ut(*o) ~ 10 6 — 10 8 cm for the typical angular 
momentum in the presupernova star found in numerical simulations (Heger, Langer & Woosley 

2 3/2 1/2 

2000). The time t satisfies t ~ 300 (^f)' 1 ( g °"^ t )" ( R ^^ > ) (^) sec, the viscous 
infall time at i? ut(£o)> where H out is the outer disk height. 

From (1), we can also write the total mechanical disk- wind luminosity as 

1 fRou t d . 9 GM NS («) 2 M otlt (t) [l- (f-) P " . 
W*) = i / T5(M dw (R,t)Vl)dR ^ V V - r> 7T\ 1 KR ° UtJ J (2) 



2j Rin dR" aw ' >-^> -~ R in (t) 

where V esc is the escape speed at the inner radius, and the last step holds when the total 
mass inside radius R satisfies M tot (R,t) < M^sR/Rui- Note that the initial total disk mass is 
~ M out (to)to, so our equations could be scaled in terms of this quantity. 



3. How p can be constrained from disk-wind theory 

For magnetically mediated outflows (e.g. Spruit 1996), the mechanical wind luminosity is 
powered by Poynting flux from large scale fields. The associated magnetic luminosity is 

L dw = — / (Ex B)2nRdR = - / V^B.B^RdR, (3) 

where E is the mean surface electric field, B z and are the surface components of the mean 
magnetic field, and is the Keplerian speed. The surface magnetic fields in (3) can be related to 
the fields in the disk by using the volume integral / V • HdV = taken over a wedge in the disk 
and converting to a surface integral. Using this, Tan & Blackman (2003) showed that if large scale 
disk fields are obtained from a helical dynamo, driven by disk turbulence and shear, the surface 
magnetic stress takes the form B Z B§ ~ Airpc^a^i 2 (H/ R), where c s is the mid-plane sound speed, 
H is the disk height and p is the mid-plane density. (The proportionality to ^pc 2 s would arise for 
any helical or non-helical dynamo, but the a ss (H/R) factor would be replaced by a ss if we used 
the averaged turbulent magnetic stress rather than the magnetic stress associated with the large 
scale fields.) Combining this magnetic stress with p = 4 J^jj V , the equation for M; n (i), the disk 
scalings (e.g. Frank et al. 2002) V r ~ a ss c s ^ ~ a ss V^^-, and plugging the result into (3) gives 



2iR,„ R 2R i (R „,/R i „)P h 
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The last equation follows from the expression for M- m (t), = {GM^s I R) 1 ^ 2 , and defining 
A = R/R m and A D = P ut/Pin- We can set (4) equal to (2), to obtain 

(1 - P )(Ag - 1) = (1 - p) (^) 2 ^ J"° X^dX = 2 ^(1 - Ag" 1 ). (5) 

Eqn. (5) shows that p is not independent of Vd w /V esc and a ss once a launch mechanism is chosen. 
The above formalism involves the relation that the mechanical disk-wind luminosity at a given 
radius is a constant fraction of the accretion luminosity, so Ref (4) implies that both L^ w and the 
accretion luminosity increase with decreasing p if all other quantities were fixed. This formalism 
breaks down when the value of p is so low that V^ w — ► c - 



4. Pulsar spin-down 

The pulsar dipole luminosity for a neutron star of angular speed f2 and radius Rns is given by 

■Bgn*(t)fl« B 

pulW 6c3 ' W 

where P* = Psin#, and 9 is the angle between the magnetic and spin axes. We have not 
considered the time evolution of sin# here. This was studied by Beskin et al. (1993) and could be 
included in future generalizations. The rotational energy evolution is determined by this dipole 
radiation and a model dependent accretion torque so that 



e = mn = -pn 4 + 2M m Rl l nn K (R m ) 



flK(Rn 



(7) 



where / is the moment of inertia of the star, (3 = B 2 sin# 2 i?^ s /6c 3 , and we use a 7 = 1 model. We 
also define the braking index and timing age as n(t) = and T(t) = ^ respectively. 

We determine the spin evolution of the star by numerically solving Eq. (7) combined with the 
equations in section 2. To do so, we also need the pulsar's magnetic field, the initial spin period, 
and the initial accretion rate in the disk. If the star were spinning down only by dipole radiation, 
then its magnetic field would be directly determined from measuring P and P. But when a disk 
torque contributes to the spin-down, we only know that the magnetic field strength must be smaller 
than the value obtained for a pure dipole. For young pulsars, we expect 10 11 ^ B <^ 10 13 G. The 
initial spin period Pi is also uncertain, but typically, 1 <J Pj <J 20 msec. The initial accretion rate is 
likewise unknown. To constrain these parameters in our model, we use observational values of the 
period, timing age and braking index for the specific pulsar under consideration. We then solve 
the coupled system of equations using a grid of values for Pj, M out (to) and R>* to obtain a best fit 
parameter range that reproduces the observed values of P, T and n to within a few percent. 



We specifically consider the Crab and Vela pulsars. In 1972, the Crab's true age was 918 yr, 
and its measured timing age was T = 1243 yr, while the braking index was n = 2.51 with period 
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P = 33.1 ms. Vela has a braking index of n = 1.4 ± 0.2 (Lyne et al. 1996), a period P = 89 
ms, and a timing age T = 11.3 kyr, for a pulsar age (determined from the associated supernova 
remnant; Aschenbach, Egger & Trumper 1995) in the range 18-30 kyr. Here we use ~ 25 kyr. 

Fig. la shows two solutions (p = 0.5 and p = 0.01) for the Crab, which reproduce the 
measured values of P, T and n. The p = 0.5 solution was obtained with = 1.5 x 10 12 G, 
Pi = 17.7 msec, and M ut(*o) = 10 29 g/s, while that with p = 0.01 has B* = 1.7 x 10 12 G, 
Pi = 17.4 msec, and M ut(*o) = 9.1 x 10 27 g/s. The bottom panel displays the pulsar and wind 
luminosities [Eq.(6) and Eq.(2)] during the pulsar lifetime. Fig. lb shows, for the same values of p, 
the solutions for Vela which reproduce its observations within the measurement uncertainties. The 
solution with p = 0.5 was obtained with B* = 7.7 x 10 11 G, Pi = 7.9 msec, and M ut(*o) = 3 x 10 30 
g/s, while that with p = 0.01 has B* = 7.6 x 10 11 G, Pi = 7.9 msec, and M ou t(*o) = 4 - 7 x 10 29 
g/s 2 Also displayed are the pulsar and wind luminosities, and the observed lower limit on the jet 
luminosity. Pavlov et al. (2003) estimated an energy injection rate ^ 8 x 10 33 g/s for the Vela jet. 
In all of our solutions, the corotation radius is less than the magnetospheric radius at all times, 
so the pulsars are always spun down by the disk. The precise value of p is not crucial to obtain a 
good solution for P, T and n; what matters for the spin-down torque is the value of M- m , and a 
different p can be compensated with a different M out to yield a similar Mj n (and hence a similar 
evolutionary track). This is why the two lines in the top three panels of Fig.l overlap. Note that 
the higher p curves for L^ w are above the lower p curves because M out is larger for the larger p 
cases, and this counteracts the dependence of Lj w on p discussed below Eq. (5). 

The particular solutions in Fig. 1 are not unique; multiple combinations of the parameters 
B, Pi and M(to) can reproduce the values of P, T and n within the observational uncertainties. 
However, we have shown solutions that also produce a disk- wind powerful enough to influence the 
jets observed in these sources, as discussed in more detail in the next section. For the solutions of 
Fig. 1, R m ~ 9 x 10 7 cm and i? out ~ 6 x 10 9 cm for Crab and R[ n ~ 10 8 cm and R out ~ 10 10 cm for 
Vela. For the Crab, the predicted optical emission from a even a face-on thin disk 3 would be below 
the optical detection (Sollerman 2003). For Vela, the observed optical flux limits require either 
(i) a disk almost edge-on to the line of sight (ii) disruption of the outer parts of the disk (maybe 
through interaction with the supernova remnant) or (iii) a geometrically thick and optically thin 
disk, like an advection dominated accretion flow (ADAF) (e.g. Esin et al. 1997). 



2 Multiplying to with M ou t(to) from our fits gives initial disk masses ~ 0.01 — 0.5M Q . The values of M ou t(io) in our 
formal solutions, like those in accretion models of gamma-ray bursts, are highly super-Eddington. Such disks would 
be dense, hot, and thick from radiation trapping (Di Matteo et al. 2002). Here the Cannizzo et al. (1990) similarity 
solutions may be inapplicable, but this does not strongly affect the present day pulsar parameter solutions. 

3 The present day disk is assumed to be optically thick, geometrically thin, and both viscous dissipation and 
irradiation by the pulsar are included (Perna et al. 2000; Perna & Hernquist 2000). 
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5. Effect of the disk-wind on the jet 

Relativistic Poynting flux pulsar winds do not easily self-collimate (e.g. Lyubarsky & Eichler 
2001); as strong outward electric forces compete with the collimating magnetic hoop stress. (This 
problem would not arise in non-relativistic outflows, as the electric field << magnetic field.) We 
suggest instead that pulsar jets may be collimated by a surrounding disk wind. 

A necessary condition for the disk wind to collimate the pulsar wind is that the ram pressure 
associated with the disk wind exceed the magnetic pressure associated with the Poynting flux 
dominated pulsar wind on scales at or less than those where the collimation is observed. Since 
the jets appear on scales >> than that of the disk, it is reasonable to make this comparison 
far outside the pulsar's light cylinder. Magnetic energy dominates the particle energy in the 
pulsar wind, so the energy outflow rate can be estimated as L pu \ ~ c(i? 2 /87r)47rr 2 , where r is a 
spherical radius. Comparing the two pressures amounts to comparing B 2 /8ir to Mdw(*)Vd w /87rr 2 . 
But B 2 r 2 = L pu i/2c, so this amounts to comparing Mdw^dw = Alw/Klw with L pu \/c, namely 
comparing the rate of change of the disk-wind momentum to that of the pulsar wind. 

Fig. 2 shows precisely this comparison for Crab and Vela. The solutions indicate that a 
collimating influence of the disk-wind on the pulsar wind is possible for the p = 0.5 case for both 
pulsars. That p = 0.5 provides higher disk-wind momenta and higher Ld w /^pul than p = 0.01 
results because the M ou t which best fits the spin parameters is larger for the large p cases. (For 
the same M out , the large p case would have a lower and momentum than the low p case, 
as discussed below Eq. (5).) Although this effect dominates, a competing effect is that Vd w (R) 
increases with decreasing p: for both pulsars, Eq. (5) for a ss = 0.1 reveals that p = 0.5 and 
p = 0.01 correspond to Vd w /V csc ~ 0.2 and Kiw/Kisc ~ 2 respectively. Collimation also requires 
that the vector sum of the pulsar wind and disk wind momenta incurs a significant vertical 
component. Even if the disk wind were spherical, when its ram pressure exceeds that of the 
pulsar wind the pulsar wind would incur some collimation. However, since the disk wind is largely 
non-relativistic, magnetic self-collimation of the disk wind is not as problematic as that of the 
pulsar wind. Self-collimation of the disk wind would favor even more collimation of the pulsar 
wind. A more detailed model of the interacting outflows is required to further develop these ideas. 
Note that Tsinganos & Bogovalov (2002) have found numerically that relativistic central outflows 
could in principle be collimated by non-relativistic disk winds, also motivating more detailed 
models. A standard limitation of such numerical calculations is the assumption of the disk field 
geometry as a given initial condition. 

At the interface of the two winds, significant dissipation can occur, and accelerated particles 
can produce the observed radiation. In addition, bending of the collimated jet could result from 
mutual non-axisymmetry of coaxial winds, or from misaligned pulsar and disk wind symmetry 
axes. Pavlov et al. (2003) suggested that a hidden wind from within the supernova remnant is 
needed to account for Vela's bent jet. 
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6. Conclusions 

Winds from fallback accretion disks below presently detectable limits may deposit enough 
momenta to collimate and bend the Crab and Vela pulsar winds and account for their observed 
jets. These disks can also account for the observed periods, braking indices, and timing ages 
of the Crab and Vela pulsars through the action of the disk torque on pulsar spin-down. Our 
results motivate further efforts to detect disks around jetted pulsars, and more detailed models of 
the nested interacting winds. Collimating a relativistic central wind by a surrounding accretion 
disk-wind may also be a relevant paradigm for active galactic nuclei and gamma-ray bursts. 

We thank L. Hernquist, K. Menou, and the referee for useful comments. 
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Fig. 1: Evolution of the period, braking index, timing age, dipole luminosity and disk- 
wind luminosity for Crab and Vela. The points mark the observed values, though 
no such constraints exist for the Crab jet power. The age for Vela is taken to be 
that of the associated supernova remnant. Boldface lines correspond to p = 0.5 
and thin lines to p = 0.01. Results are rather insensitive to the value of p. 
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